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ABSTRACT 

The declining lightcurve of the optical afterglow of gamma-ray burst GRB000301C 
showed rapid variability with one particularly bright feature at about to=3.8 days. 
This event was interpreted as gravitational microlensing by Garnavich, Loeb & Stanek 
(2000) and subsequently used to derive constraints on the structure of the GRB optical 
afterglow. In this paper, we use these structural parameters to calculate the probability 
of such a microlensing event in a realistic scenario, where all compact objects in the 
universe are associated with observable galaxies. For GRB000301C at a redshift of 
2=2.04, the a posteriori probability for a microlensing event with an amplitude of 
ATO>0.95mag (as observed) is 0.7% (2.7%) for the most plausible scenario of a flat A- 
dominated FRW universe with fim^O.B and a fraction /*=0.2 (1.0) of dark-matter in 
the form of compact objects. If we lower the magnification threshold to Am>0.10mag, 
the probabilities for microlensing events of GRB afterglows increase to 17% (57%). We 
emphasise that this low probability for a microlensing signature of almost a magnitude 
does not exclude that the observed event in the afterglow lightcurve of GRB000301C 
was caused by microlensing, especially in light of the fact that a galaxy was found 
within 2 arcsec from the GRB. In that case, however, a more robust upper limit on 
the a posteriori probability of «5% is found. It does show, however, that it will not 
be easy to create a large sample of strong GRB afterglow microlensing events for 
statistical studies of their physical conditions on micro-arcsec scales. 
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1 INTRODUCTION 

Gravitational microlensing offers a way to study the struc- 
ture of high redshift sources on micro-arcsecond scales, be- 
sides being able to constrain the mass fraction and mass 
function of compact objects in the universe. In strong grav- 
itational lens systems (i.e. systems with multiple images of 
a single background source) , the microlensing optical depth 
is of order unity (e.g. Chang & Refsdal 1979; Gott 1981; 
Young 1981). Precisely because of this high optical depth, 
these systems are perfect for resolving micro-arcsec struc- 
ture in the lensed cosmologically-distant source if it crosses 
a caustic created by the stellar-mass compact objects in the 
lens mass distribution (e.g. Chang & Refsdal 1984; Grieger, 
Kayser & Refsdal 1988; Wambsganss, Paczyriski & Schnei- 
der 1990; Wozniak et al. 2000). Moreover, because of the 
presence of multiple images, one can in principle separate 
intrinsic source fluctuations from microlensing variability. 

For example, ongoing microlensing of the lensed (opti- 
cal) images in the system Q2237-I-0305 (Huchra et al. 1985) 



has been observed ever since its discovery (e.g. Irwin et al. 
1989; Corrigan et al. 1991; 0stensen et al. 1996; Lewis et 
al. 1998; Wozniak et al. 2000). To a smaller degree and on 
longer time scales, microlensing in Q09574-561 has been de- 
tected as well (Pelt et al. 1998; see also Refsdal et al. 2000). 
The time scale of microlensing in both cases is defined by 
the relative transverse velocity between the source, lens and 
observer, which is given by the bulk velocity of the lensing 
galaxy plus the random motions of compact objects in the 
line-of-sight to the stationary quasar images, and is typically 
of order several hundred kms~^. This results in microlens- 
ing time scales of the order of months to years for solar-mass 
objects (e.g. Wambsganss 2000). 

Besides these optical sources, recently the first case 
of radio-microlensing was reported in the lens system 
B1600+434 (Koopmans et al. 2000a, 2000b; Koopmans & 
de Bruyn 2000), suggesting extremely compact relativistic 
substructure in the lensed radio source. If this microarcsec- 
scale substructure is part of a relativistic jet, the time-scale 
of microlensing variability by stellar-mass objects reduces to 
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TIME 

Figure 1 . A cartoon of how we determine the microlensing magnification as a function of time for the expanding shell source superimposed 
on the magnification patterns caused by massive compact objects. The five panels indicate five difl'erent epochs (in practise, we consider 
100 epochs). The superimposed curved line indicates the microlensing magnification as a function of time, corresponding to this particular 
example. The arrow points to a particular high magnification structure that is passed by the ring source at epoch i2. 

several weeks (Koopmans & de Bruyn 2000), allowing one 
to probe compact objects up to ^10"" Mq on time scales of 
several years, as well as the substructure of the relativistic 
jet. Indications of optical microlensing in B1600+434 with 
similar time scales have also been found (Burud et al. 2000) . 
Clearly, microlensing is a promising field of future research 
regarding the study of high-z sources at micro-arcsecond 
scales. 

In addition to testing the structure of AGNs, Loeb & 
Perna (1998) more recently proposed the use of microlens- 
ing to probe into the internal structure of GRB afterglows 
on micro-arcsec scales. Garnavich, Loeb & Stanek (2000) 
indeed interpreted an anomalous event in the lightcurves of 
the optical afterglow of GRB000301C (Masetti et al. 2000; 
Sagar et al. 2000; Berger et al. 2000; Jensen et al. 2000; 
Smette et al. 2000) as being caused by microlensing of the 
GRB afterglow. They subsequently derive constraints on its 
structure, which appear in good agreement with theoretical 
blastwave models. The inferred mass of the lensing object 
is ~0.5Mq, if its redshift is optimal for microlensing (i.e. 
about half way; Garnavich et al. 2000) . 

However, one needs to be cautious here, because no mul- 
tiple images are present - as in the case of strong gravita- 
tional lenses - to confirm that this is indeed a non-intrinsic 
event. The fact that the event occurs within only a few days 
after the burst and has an amplitude of ~1 mag, suggests 
that the burst must have occurred close to the Einstein ra- 



dius of an intervening massive compact object (Garnavich et 
al. 2000) . To have a significant probability of observing such 
a GRB microlensing event, the universe requires a surface 
density in compact objects close to the critical surface den- 
sity (e.g. Press & Gunn 1973; Blaes & Webster 1992). The 
situation is partly similar to the case of variability in single 
quasars (i.e. not multiply-imaged), which also cannot easily 
be proven to be due to microlensing (e.g. Hawkins & Tay- 
lor 1997; Hawkins 1998), although the freedom to model the 
GRB afterglow lightcurve (mostly dominated by self-similar 
expansion resulting in a power-law behavior) is significantly 
less than that for quasars. 

In this paper, we investigate the probability of mi- 
crolensing in GRB afterglows for a flat A-dominated cos- 
mological model and with the distribution of massive com- 
pact objects connected to visible galaxies, particularly focus- 
ing on the strong event seen in the GRB000301C afterglow. 
Our conclusions, however, are independent of this partic- 
ular GRB. In Section 2, we describe numerical microlens- 
ing simulations of GRB afterglow light curves, extended to 
high microlensing optical depths. In Section 3, the a pos- 
teriori probability of the observed event in GRB000301C is 
calculated. Section 4 summarizes our result and states our 
conclusions. 
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2 MICROLENSED GRB LIGHT CURVES 

To calculate the microlensing probability of GRBs, we simu- 
late microlensing magnification patterns for a range of shears 
and dimensionless surface densities in compact objects (i.e. 
7 and k; see Schneider, Ehlers & Falco 1992, Chap. 5 for the 
definitions) . A similar analysis for the microlensing variabil- 
ity in Q2237-F0305 was done by Ranch & Blandford (1991) 
and Jaroszynski et al. (1992), placing constraints on its ac- 
cretion disk models. 

As a first-order model, wo here assume that (i) galax- 
ies in the universe can be described as a singular isothermal 
sphere mass distribution (SIS; e.g. Binney & Tremaine 1987) 
for which k=^ (e.g. Kormann et al. 1994), (ii) the fraction 
of mass in compact objects /, = 1, although in Sect. 3. 3 we 
will consider the case of /*<1, and (iii) the mass spectrum 
of compact objects is narrow (Am/m<l) such that we can 
assume that all objects have nearly the same mass (for ex- 
ample the average or median value) . Magnification patterns 
for fi:=7=0.005, 0.01, 0.025, 0.05, 0.1 and 0.25 arc calcu- 
lated on a grid of 1024x1024 pixels using the ray-shooting 
algorithm of Wambsganss (1999). Each pixel has a size of 
0.1 Einstein radius. 

For the structure of the GRB, we take the results from 
Garnavich et al. (2000) , who model the GRB as an expand- 
ing ring with a width W times the ring radius R{t) . No jet 
structure is assumed in the model. The ring radius evolves 
as function of time as R{t) — i?o x t'^^ (e.g. Waxman 1997), 
where t is time in days and Rq is the ring radius on day 
one. Their best fit to the combined optical data sets sug- 
gests W^=0.16±0.02 and i?o=0.49±0.02 Einstein radius. A 
constant surface brightness inside the ring is assumed and 
no emission comes from either inside or outside of the ring 
structure. 

Although the redshift and mass of the possible lensing 

object is unknown, this is of no consequence to our calcula- 
tions, because we can express the properties of lensing ob- 
jects in dimensionless units of critical surface density, shear 
and Einstein radius. For intermediate redshifts the mass of 
the possible lensing object of the GRB000301C afterglow 
would be --^0.5 Mq (Garnavich et al. 2000), which agrees 
well with Galactic (e.g. Alcock et al. 2000) as well as ex- 
tragalactic constraints on halo mass objects (Refsdal et al. 
2000; Wambsganss et al. 2000). 

To obtain simulated GRB microlensing light curves, we 
convolve the magnification patterns with the time-variable 
GRB source structure for each epoch t. We then store the 
'convolved' magnifications from lO'' random grid points. We 
repeat this procedure for 100 epochs from day 0.5 to day 
35 (approximately the period during which the afterglow 
lightcurvo of GRB000301C was sampled; Garnavich et al. 
2000), each step evolving the GRB afterglow structure and 
taking the magnifications from the same 10* grid points (the 
procedure is illustrated in Figure 1). The sampling epochs 
for our simulated microlensing lightcurves are logarithmi- 
cally spaced, to avoid undersampling during the initial GRB 
afterglow phase, where one might expect the highest magni- 
fications. The light curves are normalized to show only the 
microlensing magnification of the GRB afterglow (as in Mao 
& Loeb 2000). In Figure 2, we show the probability density 
distributions of the peak magnification of the normalized 
GRB microlensing light curves, from which it is immediately 



K=7 


P(Am>0.95™) 


P(Am>0.10™) 


0.005 


0.0043±0.0007 


0.181±0.004 


0.010 


0.0063±0.0008 


0.224±0.005 


0.025 


0.0150±0.001 


0.441±0.006 


0.050 


0.0034±0.002 


0.617±0.008 


0.100 


0.05.30±0.002 


0.729±0.009 


0.250 


0.0530±0.002 


0.781±0.009 



Table 1. The a posteriori probability that GRB000301C shows 
a microlensing event with magnitude greater than approx. 0.10 
and 0.95 mag, respectively. The errors indicate the Poisson error, 
due to the finite number (i.e. 10*) of simulated light curves. The 

equal probability for «;=0.1 and 0.25 in the second column is 
coincidental, but illustrates the rapid break from Pock (see text) 
in the probability of strong events for high optical depth regimes. 



clear that strong fsl mag events are relatively rare. We will 
discuss this in more detail in the next section. 



3 PROBABILITIES 

In this section we consider the probability of observing the 
event seen in GRB000301C, within the microlensing hypoth- 
esis. First, we calculate the probability of the event as func- 
tion of the dimensionless surface density in compact objects. 
Second, we determine the probability of actually observing 
GRB000301C through a dimensionless surface density, k, as 
function of a cosmological model and distribution of com- 
pact objects. Third, we combine these probabilities to arrive 
at an a posteriori probability for the feature in the afterglow 
light curve of GRB000301C being a microlensing event, as 
suggested by Garnavich et al. (2000). 



3.1 Event- Amplitude Probability 

The microlensing event seen in GRB000301C has a maxi- 
mum amplitude of 0.95 mag or a magnification /io=2.4 (see 
Figure 2 in Garnavich et al. 2000). To calculate the proba- 
bility of at least such a strong event, we need to measure the 
fraction of simulated light curves (Sect. 2) that show an event 
stronger than /Xo (see Figure 2). These numbers are listed 
in Table 1. We find that even in the higher optical depth 
regimes the probability does not significantly exceed ~5%, 
but levels off in the high optical depth regime, because of av- 
eraging of the GRB afterglow over multiple caustics. At low 
surface densities, one expects the probability of microlensing 
(r^i) to asymptotically behave as r^iocK. 

Let us now make a simple analytical fit through these 
probabilities, using the following function 



PcK 



(1) 



This function is by no means the correct functional form, 
that might be expected from a detailed theoretical analysis, 
but it has the correct asymptotic behaviour for k— >0. We 
find that Pc«0.82 and Kc~0.072 fit the results in Table 1 
best. For the purposes here, this empirical function provides 
a good enough representation of the probability of the ob- 
served event in the region O < /t < 0.25. The value /t=0.25, 
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is the highest through which any source can be seen with- 
out being multiply-imaged in case of a SIS mass distribution 
(e.g. Kormann et al. 1994). In Table 1, we also list the re- 
sults for events with amplitudes greater than 0.1 mag (i.e. 
^o = l-l). In that case, we find Pc~22.7 and We 
defer a discussion of these events to Sect. 3. 3. 

For the high-magnification events (>0.95mag), Table 1 
shows that P^k. Because the probability of multiple imag- 
ing by compact objects is equal to k (in the case of k<^1), 
this result implies that the GRB uccds to lie very close 
to the Einstein radius of the compact object as found 
for GRB000301C (Garnavich et al. 2000). Conversely, one 
could conclude: The fact that the microlonsing event in 
GRB0QQ3Q1C requires the GRB to lie close to the Einstein 
radius, implies a probability of the event close to the av- 
erage dimensionless surface density <k> of the universe in 
compact objects. 



3.2 Surface Density Probability 

To obtain the final probability of this event being observed, 
we have to multiply eqn.(l) with the probability that the 
GRB is seen through a region of sky with a dimensionless 
surface density between k and K+dn. 

For a constant comoving density of SIS lens galajc- 
ics, which distribution follows the Schcchtcr luminosity (i.e. 
mass) function, we find from Turner, Ostriker & Gott (1984), 
Turner (1990) and Fukugita & Turner (1991) that the opti- 
cal depth for multiple imaging is 



-rGL(2:s) 



VHo 



a H hi 

7 



G{z.), (2) 



where and a*j are the local density and central velocity 
dispersion of L* galaxies, respectively. The parameters in 
the F-function describe the Schechter luminosity function 
(see Fukugita &i Turner (1991) for a detailed description of 
this equation). Furthermore 



dw 



(3) 



for a fiat FRW universe with fim -|- f2A=l and a GRB red- 
shift ^s. For definiteness, we assume here that f2m=0.3 and 
f2A=0.7, which seem to agree best with most recent CMB, 
SNe la and cluster-abundance observations. 

The optical depth tgi,{zb) indicates the fraction of sky 
covered by regions in which sources are multiply imaged (i.e. 
regions inside the Einstein ring). It is easy to show that the 
brightest imago, in case of multiple imaging of the source, 
always lies in the region k=0.25-0.50 (between 1-2 Ein- 
stein radii from the lens centre) for a SIS mass distribution 
{k—1/2x, with X being the radius in Einstein radii). This 
region projects one-to-one back onto a disk of one Einstein 
radius in the source plane (see Schneider et al. 1992, Chap- 
ter 8). Hence the probability of observing a GRB through 
a patch of sky with dimensionless surface density >k (but 
smaller than «;=0.5) is given by 



P(> ^) = rGUz.)x (^-l)' 



(4) 



in which case the region /t=0.25-0.50 exactly has the proba- 
bility tgi^{zb), as required. We do not take the magnification 



bias into account here, because the majority of GRB will 
be seen through regions with n<^l. In Sect. 4, we will come 
back to the magnification bias in more detail. Finally, we 
need to normalize P{> k) to unity, which requires us to put 
a lower limit on the allowed k. We find fi;i=|(l + I/^tgl)^'' 
for which P{> Ki)=1. For surface densities smaller than 
Ki, galaxies will start to significantly overlap. In that case 
eqn.(4) will clearly break down and Ki is therefore a natural 
limit on the probability distribution of k. 



3.3 Probability of the GRB000301C Event 

The overall probability that a GRB with the properties of 
GRB000301C shows a microlensirig event with an amplitude 
of RiO.95 mag (i.e. /io=2.4) then becomes 



r 

JO.5 



dP(> k) 
dn 



X PA{K,\Ho)dK. 



(5) 



If we now use the values no=0.61xlO~^ft.^Mpc~^, a=— 1, 

7=4 and aj^=225 km/s for the Schechter luminosity function 
describing the population of elliptical lens galaxies, which 
presumably dominates the lousing cross-section and mass 
in the universe (see also Kochanek 1996; Falco et al. 1998), 
we find 



rGL(2s) w 9 • IQ-* X G(0s). 



(6) 



Evaluating eqn.(6) for the redshift of the burst 2:grb=2.04 

(Jensen ot al. 2000; see also Smotto et al. 2000), we find 
the probability that GRB00301C could have been multiply 
imaged to be 



tgl/grb 



2.2 X 10" 



(7) 



for Oni=0.3 in a fiat FRW universe. Hence the probability 
of seeing the observed microlensing event becomes 



Pfi ^ tgl/grb 



(2 K - 1) Pc 



0.5 2«2 V^l -h [K/Kcf 



(8) 



For GRB000301C at ^^=2.04, we furthermore find ki«0.022 
(Sect. 3.2). Evaluating eqn.(8) for Pc^O.82 and Kc~0.072 
(Sect. 3.1), we finally find that the a posteriori probability 
of the event seen in GRB000301C is Pe(/« = 1) w 0.027. 

However, this calculation assumes that all mass is in 
the form of compact objects. Recent results from the MA- 
CHO collaboration indicate that /, =0.08 0.50 with 95% 
confidence (Alcock et al. 2000). The most likely value is 
/tRiO.2. In that case, we have to modify eqn.(l) such that 
K^/, K. In other words, the probability of an event goes 
down by a factor ~/*. This scaling of the numerical simu- 
lation is allowed, because for «;=7<Cl, the infiuencc of the 
shear on lensing properties is small (see also Mao & Loeb 
2000). Similarly, changing a fraction (1 — /») of the surface 
density in compact objects to a smooth mass distribution is 
nearly similar to completely removing this fraction for fv<Cl. 
Thus for a more realistic scenario, the probability reduces 
to Pe(/, = 0.2)w0.007. 

For microlensing events with an amplitude Am > 
O.lmag (much weaker events than the one seen in the 
GRB0003001C afterglow), we find Pe«0.56 and 0.17, re- 
spectively, for /,=1.0 and 0.2. One might therefore expect 
many high-z (i.e. z^2) GRB afterglows to show evidence for 
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Figure 2. The probability distribution function of microlensing magnifications as function of k=7. From upper left to lower right, 
K=7=0.005, 0.01, 0.025, 0.05, 0.1 and 0.25, respectively. The dashed line indicates a peak magnification of 0.95 mag or /^o=2.4. 



microlensing at a low level, assuming their physical prop- 
erties are similar to those inferred from GRB000301C by 
Garnavich et al. (2000). 



4 DISCUSSION 

If the event seen in the optical lightcurve of the 
GRB000301C afterglow is caused by microlensing, the op- 
portunities to study the structure and evolution of GRB 
afterglows on micro-arcsec scale are potentially very excit- 
ing (Loeb & Perna 1998; Garnavich et al. 2000; Mao & Loeb 
2000). 

However, in this paper we have shown that the a pos- 
teriori probability of this particular event is actually small, 
i.e. 0.7-2.7% for a fraction /,=0. 2-1.0, respectively, of dark- 
matter in the form of compact objects. The main assump- 
tions in this calculation are: (i) a constant comoving density 
population of galaxies which follow a Schechter luminosity 
(i.e. mass) function, (ii) a flat A-dominated FRW universe 
with S7ni=0.3, (iii) all matter in the universe traces galaxies, 
which can be described as singular isothermal spheres, (iv) 
the mass spectrum of compact objects is 'narrow', in which 
case they can be parameterised by a single value for their 
mass, (v) a typical GRB optical afterglow has similar prop- 
erties as GRB000301C and has a redshift of z^2, and (vi) 
there is no significant magnification bias. The magnification 
bias, however, is unlikely to be a problem. Even in case this 
bias increases the number of observed GRBs by a factor of 
10 for high microlensing optical depths (k<;0.25), Table 1 



and eqn.(7) show that the overall probabilities are increased 
by only 10xP(Am > 0.95'^\k ^ 0.25) x rGL/GRB~0.1%. 
A detailed calculation should consistently take into account 
both the magnification distribution of galaxies and that of 
compact objects (e.g. Pei 1993a, 1993b). In light of the very 
uncertain redshift distribution and luminosity function of 
GRB afterglows this is, however, not yet warranted. 

Moreover, even in the case that the whole sky has a sur- 
face mass density of n,^0.25 (normalised by the "critical" 
surface mass density for lensing, e.g. Schneider et al. 1992) 
our simulations (see Table 1) indicate that the probability 
of such an event does not significantly exceed ~5%. Clearly, 
the latter situation is unrealistic, but it does place a very 
robust upper limit on the probability of this event (see also 
Press & Gunn 1973 and Blaes & Webster 1992). 

The plausibility to actually observe such a fraction of 
GRB afterglow lightcurves with a microlensing event of 
about one magnitude also depends on the typical mass of 
compact objects. However, even if the mass goes up or down 
by a factor of ten from the assumed 0.5 M0, to first approx- 
imation that would just mean a "broadening" in time of 
the microlensing event by about a factor of three (since the 
length/time scale is proportional to the square root of the 
mass). However, if most compact objects have masses much 
smaller than that inferred for GRB000301C, the fraction of 
GRB afterglows with microlensing signatures will obviously 
decrease. In that case the effect of the "shrunk" caustics will 
average out over the much larger angular size of the emis- 
sion region of the GRB afterglow. This case, however, seems 
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unlikely in light of the stringent lower limit (^lO'^M©) 
placed on low-mass compact objects in 0957+561 (Schmidt 
& Wambsganss 1998; Refsdal et al. 2000; Wambsganss et al. 
2000). The fraction of microlensed GRB afterglows would 
also decrease if the typical redshift of GRBs is significantly 
lower than z = 2. A somewhat lower typical redshift is sup- 
ported by the average redshift < z >grb~ 1-3, with an rms 
spread of «1, that we find from Bloom, Kulkarni, & Djorgov- 
ski (2000) . In case of a broad mass spectrum of the compact 
objects or a large spread in the physical properties of the 
GRB optical afterglows, one has to convolve the microlens- 
ing probabilities with these properties. This, however, is not 
expected to change the results significantly. 

We emphasise that our results do not exclude that 
the event seen in the lightcurve of the optical afterglow of 
GRB000301C is in fact due to gravitational microlensing, 
especially in light of the fact that a galaxy was found at 
2 arcsec from the GRB (Garnavich et al. 2000). Each event 
should be treated on its own merit (and it is a posteriori 
anyway). In any case, in the sample of ~20 known GRBs 
with observed optical afterglows (e.g. Bloom et al. 2000), 
the probability of one such event would be around 14%. 
But our results do predict that it is unlikely to find many 
strongly microlensed GRB afterglow soon, and it emphasises 
the difficulties one will encounter in creating a large sample 
of these strong (i.e. Am ^ 1 mag) events, with the aim to 
statistically study the physical properties of GRB afterglows 
on micro-arcsec scales. For example, to get 10 additional mi- 
crolensing events of comparable strength, one would require 
~1500 GRB afterglow optical lightcurves. Even for the Swift 
satellite (Parsons et al. 1999; see also Mao & Loeb 2000) this 
provides a challenging task. 

On the bright sight, if the event in GRB000301C was 
due to microlensing, we can expect a significant fraction 
of high-z GRB afterglow light curves to show microlensing 
events stronger than 0.1 mag (see also Mao & Loeb 2000), 
although these will be hard to distinguish from intrinsic vari- 
ations and will also not really be able probe the GRB af- 
terglow structure on micro-arcsec scales. We conclude, that 
the best cases to unambiguously study micro-arcsec struc- 
ture in GRB afterglows will probably be those GRBs that 
are multiply imaged (about 1-2% expected, see, e.g. Holz 
et al. 1999; Romberg et al. 1999; Mao 1992, 1993; Marani 
et al. 1999; McBreen et al. 1993; Narayan & Wallington 
1992; Nowak & Grossman 1993; Paczyhski 1987; Wambs- 
ganss 1993; Williams & Wijers 1997). Those (soon to be 
discovered!) will almost certainly (see Table 1) show mi- 
crolensing events ^0.1 mag, which can easily be separated 
from intrinsic variations by comparison with the other lensed 
GRB images, after correcting for the respective time delays. 
A smaller fraction (~5%) will also show events ^1 mag, al- 
though this might increase due to the magnification bias. So 
the prospects of doing science with (micro-) lensed gamma 
ray bursts and afterglows remain promising. 
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